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Collisional avalanches are identified to be responsible for an 8-fold increase of the initial loss rate 
of a large 87 Rb condensate. We show that the collisional opacity of an ultra-cold gas exhibits a 
critical value. When exceeded, losses due to inelastic collisions are substantially enhanced. Under 
these circumstances, reaching the hydrodynamic regime in conventional BEC experiments is highly 
questionable. 
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One of the current goals in the field of Bose-Einstein 
condensation (BEC) is the production of a condensate in 
the collisionally opaque or hydrodynamic regime, where 
the mean free path of an atom is much less than the size 
of the sample. This would offer the opportunity to study 
striking phenomena like quantum depletion or dynami- 
cal local thermal equilibrium. In this context, one pos- 
sible approach is to increase the interaction among the 
atoms by means of Feshbach resonances pj . It has been 
observed, however, that in their vicinity the large cross- 
section for elastic collisions is accompanied by a dramatic 
increase of atom losses 0, . Hence, it seems advanta- 
geous to follow a different route by producing large and 
dense condensates. 

In this letter we conclude that the collisionally opaque 
regime can hardly be reached in alkali BEC experiments. 
We identify an intrinsic decay process that severely limits 
the average column density (nl) of condensates at values 
achieved in present BEC experiments. It is based on 
collisional avalanches that are triggered by inelastic col- 
lisions between condensate atoms. A considerable part 
of the energy released in these initiatory collisions is dis- 
tributed among trapped atoms resulting in a dramatic 
enhancement of the total loss from the condensate. In 
analogy to the critical mass needed for a nuclear explo- 
sion, we define a critical value of the collisional opacity 
(nl)a s , with a s = &ira 2 the s-wave cross section for like 
atoms and a the scattering length. The critical opacity 
equals 0.693, corresponding to a collision probability of 
0.5. Related scenarios have been discussed in Refs. 
but were assumed to play a minor role in the experimen- 
tally relevant region. However, we present strong exper- 
imental evidence that the anomalous decay of our 87 Rb 
condensate is caused by collisional avalanches. This is 
supported by the good agreement of a simple model with 
the data. 

The crucial point for the occurrence of an avalanche is 
whether the products of a one-, two- or three-body decay 
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FIG. 1: Sketch of a collisional avalanche in a homogeneous 
condensate with a radius equal to twice the mean free path 
(nas)^ 1 , illustrating the enhancement of the loss rate. 



process have a substantial probability 

p{E) = l-exp[-(nl)a(E)}, 



(1) 



of undergoing secondary collisions before leaving the trap 
[6j, with a(E) the total cross section at kinetic energy E. 
The collision probability varies significantly with tem- 
perature and is usually highest in the s-wave scattering 
regime. Here, the differential cross section is isotropic 
in the center-of-mass system; in the lab system the two 
atoms fly apart at an angle of tt/2 on average. The en- 
ergy of the projectile is on average equally distributed 
among the two colliding atoms. This implies that each 
collision results in two new atoms that both can continue 
their collisional havoc in the trap until they leave the 
condensate (Fig. QJ. If the probability for collisions is 
higher than 0.5, the average number of colliding atoms 
increases with every step of the collisional chain which 
now becomes self-sustaining. 

To calculate the total loss from the condensate, we 
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start from the well known loss rates Ni — —KiN^n 1 ^ 1 ), 
with i — 1,2,3 associated with one-, two- and three- 
body decay processes with rate constant i-Q, respectively. 
Here, N is the number of atoms in the gas with the den- 
sity distribution n(r). Depending on the energy of the 
decay products, typically a few or even no further colli- 
sions are needed to generate an atom with an energy Ei tS 
whose next collision would be in the s-wave regime. The 
probability for this collisional chain is pi,i ■ pi^. • • • = Pi 
with pi^ n — p(Ei tn ). During this process, on average 
§i atoms are lost from the condensate without under- 
going secondary collisions. Next, an atom with energy 
E is induces an avalanche with a collision probability 
p s = 1 — cxp[— (nl)a s ] that now is independent of en- 
ergy. Consequently, atoms with an energy E is /2 k in the 
fe** 1 step of an avalanche are generated with a probability 
of p~ip k s - Since every collision now results in two projec- 
tiles in the next step, the degeneracy of step k is 2 k . The 
rate at which atoms are lost from this avalanche step is 
PiNi2 k p k (l —p s ) and the total loss rate from the conden- 
sate becomes 
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where the sum extends over all relevant avalanche steps. 

To determine the cutoff fc|™ ax , note that avalanche- 
enhanced losses can occur up to a step with an energy 
on the order of the chemical potential. However, when 
the energy falls short of the trap depth, E tra p, the atoms 
lost from the cold sample are still trapped. They will re- 
peatedly penetrate the cloud and thus give rise to heat- 
ing. This will either be compensated by an evaporation 
of atoms from the trap or will reduce the condensed 
fraction by increasing the temperature. Both possibil- 
ities are not described by Eq. (J2J). Therefore, we use 
fcrnax <- i g 2 (^E is / E trap ) as a cutoff and, hence, account 
only for immediate trap losses. 

The additional heat-induced depletion caused by 
trapped avalanche atoms can easily be estimated for 
the case that the temperature is fixed by the trap 
depth. Each atom participating in step (kl nax + 1) of 
the avalanche will finally dump about the energy E tra p/2 
into the system. Since any evaporated atom takes the 
energy E trap with it, about half as many atoms as are 
produced in the step (k r ( lax + 1) of the avalanche have to 
be evaporated to keep the temperature constant. Hence, 
the evaporation rate is 



K,heat - (1/2) PiNi 2 
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Equation J5J predicts substantially enhanced losses as 
soon as the critical opacity is exceeded. However, for a 
given k™ ax there is a second critical value of the opacity 
above which the loss rate Ni, ava i decreases again. Now, 



with increasing opacity the limited trap depth continu- 
ously looses its shielding effect against the products of in- 
elastic collisions, since most avalanches generate trapped 
particles. In the collisional regime with p s ~ 1, the en- 
ergy released in an inelastic process will be entirely dis- 
sipated in the system. This results in an explosion-like 
particle loss according to Eq. 

To apply our model, the column density must be eval- 
uated according to 
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[n(f)/N] ■ J [n(r + R)/4ttR 2 } d 3 Rd 3 r (4) 
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dx(l+x 2 )- 1 (l+e 2 x 2 )- 1 ^ 2 (5) 



where the second line is the result for the case of a har- 
monic potential with cylindrical symmetry. Here, e = 
<jj\\/<jj±_ is the ratio of the trap frequencies and n p is the 
peak density of the cold sample. For the parabolic den- 
sity distribution of the condensate, W± is the half radial 
width and c = 5/12. Due to the scaling n p oc TV 2 / 5 and 
W± cx N x / b the column density of a condensate scales as 
iV 3 / 5 , so that the effect of multiple collisions is quite per- 
sistent. For a Gaussian distribution we find W± = a±, 
c = \JtF[% and a scaling according to (nl) oc N. In a 
harmonic potential, the ideal Bose distribution can be 
represented as a sum of Gaussian distributions and the 
latter result can thus be used to evaluate the column 
density close to degeneracy. For a Bose distribution the 
opacity scales disproportionate to N, resulting in a faster 
decline of the avalanche enhancement. 

The next step is to identify the energies of the ini- 
tial decay products. For a background gas collision, E\ \ 
depends on the mass of the impinging particle that is as- 
sumed to be Rb in our system. In the case of spin relax- 
ation, £2,1 equals either the Zeeman energy or the hyper- 
fine splitting energy. For three-body recombination £3.1 
has to be derived from the binding energy of the most 
weakly bound level in which the dimer is predominantly 
formed. Clearly, the molecule is likely to be deactivated 
in a subsequent inelastic collision with a condensate atom 
0>I3- Deactivating collisions will be a serious problem in 
highly opaque clouds where atoms with higher energies 
still have high collision probabilities. In our experiment, 
however, the collision probability is significantly smaller 
at typical deactivation energies of 0.1 K than at the bind- 
ing energy of the molecules in the last bound level. In 
our analysis we therefore do not account for avalanches 
triggered by deactivating collisions. The values of all pa- 
rameters used to calculate the effective losses are listed in 
TablelJ Note that in order to account for the avalanches 
triggered by the two-body decay, the partial rates associ- 
ated with the various exit channels are needed since they 
correspond to different energies released in the process. 

Finally, the presence of a diffuse atom cloud in the 
trapping volume can cause additional losses (see e.g. Q). 
In a steep magnetic trap with a depth of a few mK, such 
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TABLE I: Rate constants for the initiatory processes and en- 
ergies of the subsequent collisions that are necessary to gen- 
erate the first avalanche atom with energy E s . 
i type rate constant AEi , . . . , AE n = E s 





K t 




[(3/2)**] 


1 background 


1/(39 s 




4K, 100, 5, 0.5mK 


2 Zeeman b 


1.4 x 10" 18 


cm 3 /s 


0.022 mK 


2 x Zeeman c 


3.7 x KT 17 


cm 3 /s 


0.045 mK 


hyperfine d 


2.2 x 10~ 16 


cm 3 /s 


109, 5, 0.5 mK 


2xhyperfine e 


1.3 x 10 -16 


cm 3 /s 


219, 8, 0.5 mK 



3 recombination-^ 1.8 x 10" 29 cm 6 /s 0.54 mK 9 

a IH; ^^mill; 9 |H; 



an "Oort" cloud is mainly a consequence of incomplete 
evaporation at high magnetic fields |9j or low radio fre- 
quency (rf) power. In our experiment, the temperature of 
the diffuse cloud will probably be on the order of 400 /iK, 
corresponding to the measured initial temperature of the 
magnetically trapped cloud. Even in a rf-shicldcd trap 
these atoms will penetrate the condensate giving rise to 
an additional decay rate according to 1/r = n OOIt a s v OOIt , 
with n oort and iJ 00 rt the density and the thermal velocity 
of the penetrating atoms, respectively. Collisions with 
Oort atoms will also trigger avalanches, because the col- 
lision energy is close to the s-wave regime. 

To compare our data with the predictions of the model, 
the differential and the total scattering cross sections are 
needed. Above a kinetic energy E/ftB of 60 mK we calcu- 
late the energy transfer by collisions using a model func- 
tion for the small-angle differential cross section [lfj] . For 
collisions below 60 mK we use the numerical results from 
a full quantum treatment ^lj. For 87 Rb in the |2, 2} 
state, the large contribution of a d-wave scattering reso- 
nance to the total cross section leads to a(E) ~ 4 x a s 
at an energy of E/ks = 560 /iK in the lab system. This 
almost exactly coincides with the energy transferred to 
the third atom in a recombination event (tableP). Hence, 
a secondary collision of this atom will occur with a prob- 
ability of 0.99 already when the probability for s-wave 
collisions among condensate atoms is 0.7. For kinetic 
energies E/ks < 1-5 mK, the total cross section obeys 
<j(E) > <7 S . For simplicity, we use a s for calculating 
avalanches in this energy range. Our model therefore 
yields a lower bound for the total loss. 

The apparatus used to study the condensates has been 
described previously 0, The experiment is per- 

formed with 87 Rb atoms in the |2, 2} state. A Ioffc- 
Pritchard magnetic trap with a bias field of 2 G and oscil- 
lation frequencies of uj±/2it = 227 Hz and uj^/2tt = 24.5 
Hz is used. The atoms are cooled by rf evaporation and 
then held in the trap for a variable time interval. During 
the storage time, the trap depth is set to E tra p/kB = 4.4 
/iK by means of the rf shield. From the width of the den- 
sity distribution after expansion the atom number Nc in 
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FIG. 2: Decay of the condensate and the thermal cloud. The 
horizontal line corresponds to the critical opacity. For com- 
parison, the calculated decay due to the initial one-, two- 
and three-body loss rates without (dotted line) and with 
avalanche-enhancement (dashed line) are shown. The full line 
includes the effect of an Oort cloud with avalanche enhance- 
ment. 

the condensate is determined. At minimum storage time, 
we find Nc = 1-1 x 10 6 atoms and n p — 6.4 x 10 14 cm' 3 
and no discernible non-condensed fraction. 

The decay curve of the condensate is shown in Fig. 
revealing that about half the initial number of atoms is 
lost within the first 100 ms. The dotted line shows the 
theoretical prediction assuming that losses occur solely 
due to background gas collisions, spin relaxation and re- 
combination (Table[IJ|. The observed loss is 8 times faster 
than predicted. Moreover, the additional decay is clearly 
non-exponential and can therefore not result from pri- 
mary collisions with Oort atoms. Hence, multiple colli- 
sions have to be taken into account. 

Indeed, with (nl)a s = 1.4 the critical opacity is consid- 
erably exceeded. To the best of our knowledge, the cor- 
responding s-wave collision probability of p s = 0.76 has 
not been reached in published work on Rb condensates in 
the off-resonant scattering regime. This explains why our 
observations differ from those made in other experiments 
0003- The dashed line displayed in Fig. |2|has been 
obtained by numerically integrating the rate equation 
N = Y^i=i(^i,aval + NiMeat) that describes avalanche- 
enhanced losses according to Eqs. J5J) and ©, without 
any adjustable parameter and neglecting the contribution 
of an Oort cloud. We find good agreement between the- 
ory and experiment within the first 200 ms, showing that 
collisional avalanches triggered by recombination events 
are responsible for the fast initial decay. 

To investigate the role of an Oort cloud, we have 
performed a similar experiment with an atom cloud at 
a lower density. Figure shows the decay of a non- 
condensed cloud with 1 x 10 7 atoms at a temperature of 
1 /iK and a peak density of 3.5 x 10 14 cm -3 . The num- 
ber of atoms is determined from the total absorption of 
near-resonant laser light. The trap depth is limited to 
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10 /iK, according to the higher temperature of the sam- 
ple. Again, the decay is non-exponential and initially 
about two times faster than predicted by the primary 
loss-rates (dotted line). At an opacity of 0.9, obtained by 
assuming an ideal Bose distribution, we already expect 
a weak avalanche-enhancement. This allows us to test 
our model in a different regime since in a thermal cloud 
avalanches are less persistent than in a condensate. In ad- 
dition, the intrinsic two- and three-body decay rates will 
die out during the observation time whereas the effect 
of an Oort cloud as a one-body decay will persist. The 
solid line in Fig. [2]is the prediction of our model where we 
have included an avalanche-enhanced decay rate caused 
by an Oort cloud. Good agreement with the data is ob- 
tained for 1/r = 1/7. 8s, corresponding to n oor t = 5 x 10 8 
cm~ 3 at 400 fiK. Such a density is produced by only a 
few times 10 5 atoms and appears realistic in view of the 
more than 10 9 atoms that were loaded into the magnetic 
trap. It is also consistent with the fact that we have no 
direct experimental evidence for an Oort cloud and that 
the initial decay is correctly predicted by the model even 
if the contribution of the Oort cloud is neglected (dashed 
line) . 

We can now calculate the extra loss rate of the conden- 
sate due to an Oort cloud. Since the two experiments de- 
scribed above are performed under identical conditions, 
the density of the Oort cloud is essentially unchanged in 
the two measurements. As can be seen from the solid 
line in Fig. [21 the small extra contribution from the 
Oort cloud does not change the predicted initial decay 
but slightly improves the agreement between the model 
and the data for longer times. The small remaining dis- 
crepancy can be the result of an additional decay not 
accounted for in our model. In particular, avalanches 
will seriously perturb the equilibrium of the condensate 
by inducing local fluctuations of the mean-field energy 
Since the damping rate of excitations can be small 
compared to the elastic collision rate, we expect that this 
process introduces a second time scale to the decay that 
depends on the history of the condensate. 

The simultaneous agreement of our model with the two 
complementary data sets strongly supports the evidence 
for the occurrence of collisional avalanches in our exper- 
iments. Our analysis reveals that the density of a cold 
gas is severely limited as soon as the s-wave collisional 
opacity exceeds the critical value of 0.693. It is impor- 
tant to point out that the anomalous initial decay of the 
condensate is attributed to collisional avalanches almost 
exclusively triggered by the intrinsic process of recom- 
bination and that no free parameters are introduced in 
the model. We have no evidence for the contribution of 
an Oort cloud to the fast initial decay observed in our 
experiments. 

We conclude that it will be hard to enter the collisional 



regime in alkali BEC systems. For 87 Rb in the \2, 2) state 
the prospects are even worse due to the large collision 
cross section of the recombination products. Hydrody- 
namic conditions might be reached in the longitudinal di- 
rection in an extremely prolate geometry, as can be seen 
from Eq. JSJ. Avalanches triggered by recombination 
events can be suppressed in the vicinity of Feshbach res- 
onances, where the recombination energy becomes very 
small. However, collisional deactivation of the highly ex- 
cited molecules can still produce avalanches Q and might 
contribute to the fast decay reported in Ref. 0. This 
offers a new application for a condensate of, e.g., ground- 
state helium atoms, where recombination is not possible. 
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the scattering cross sections and the spin relaxation 
rates. H.C.W.B. gratefully acknowledges the hospitality 
of JILA, NIST and the University of Colorado, Boulder. 
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